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Abstract 


A kinetic expression for the chemical reaction of lithium metal with thionyl chloride is presented that is consistent with calorimetric measurements 
of the heat generation from a thionyl chloride cell. The kinetics expression is incorporated into a well-established electrochemical model for the 
discharge behavior, and then used to estimate the life of the battery under an intermittent discharge so as to assess the importance of lithium 
corrosion. The model predicts that, under the conditions examined, there is no danger of depleting the lithium anode and so introducing a safety 


hazard. 
Published by Elsevier B.V. 
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1. Introduction 


Self-discharge in Li/SOCl, batteries is largely associated 
with the reaction of the lithium anode with the thionyl chlo- 
ride electrolyte (see reaction (1) below). Estimation of the self- 
discharge of Li/SOCl, batteries is an important problem in the 
design of safe cells and batteries. If self-discharge leads to a 
complete loss of lithium anode in one cell in a string of cells, 
then that cell can go into reversal when the string is discharged 
and possibly explode. Thus, it is important to estimate the self- 
discharge rate of the lithium anode. Long-lived (~10 years) 
thionyl chloride batteries are of interest for applications, such 
as utility meters and other applications. Over such long periods 
of time, the self-discharge reaction should become more sig- 
nificant. Despite its importance, the issue of self-discharge has 
received little attention compared to the discharge behavior. This 
paper presents an approach for estimating the self-discharge of 
the lithium anode under various pulse discharge scenarios. 

The design and general of behavior of Li/SOClp batteries is 
described by Linden and Reddy [1]. The desired electrochemical 
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reactions are, 


Anode: Li > Lit+e7 


Cathode: 4Lit +4e7 + 2SOCl, > 4LiCl + S + SO, 


A key design requirement is that the cathode must be porous 
enough to accommodate the solid LiCl formed on discharge. 
This requirement complicates the design because the cathode 
can expand during discharge. Although several groups have 
reported models to simulate discharge of Li/SOCl cells [2-5], 
the problem of cathode expansion has not been addressed. Others 
have considered the thermal [6-8] behavior of Li/SOCl batter- 
ies. However, the problem of self-discharge has received little 
attention. 

The self-discharge reaction, 


2SOCl2 + 4Li > 4LiCl + S + Ss (1) 


is known to form a lithium chloride film at the lithium anode. 
Roth [9] characterized the self-discharge of Li/SOCl2 D-size 
cells using a CSC4400 model microcalorimetry that was cal- 
ibrated with heaters of similar size and shape to D-cells and 
corrected the results for background heat leaks. Roth posited 
that the heat release was attributed to reaction (1) and used 
the enthalpy change of that reaction (AH =—3.65E5 J mol~!) 
to estimate the lithium loss over time. Roth found that the rate 
of reaction (1) depended on how the cell was used. Immedi- 
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Table 1 
Summary of model equations and parameters 
Name Equation 
(a) Equations 
S m HE A 
Kirchoff’s law — = Fjy 
ox : 
3 an, 
Solid-phase Ohm’s law Zë SE 
dr Gett 
d j 2RT Ve d 
Liquid-phase Ohm’s law $2 Zara ( 2 EO = : 
dr Kë cF 21—cV)/ dr 
an Sei dec a. eDac SIS À ey, Wa i2, x Oty 
alt mass balance = t+ CV] + Jy i 
ax Lt Ox i PAT ia eae E a? 
; an ba y , Qa 1 F c =e 1 F , 3 —4641 
Negative electrode kinetics iy = 10,1,ref | €XP nı exp nı > nı = di $2 — Ui sets io,1,ref = 1.157 x 10° exp 
RT Cref RT 


Qa 2 F 


Positve electrode kinetic Jy = 10,2,retAs fe ( 


2 exp n2 > SLEP] zg 
RT Cret / A Onst RT 


Volume balance 2 = Viicijy 
Darcy’s law v= [a OH) (Yo = sa E 
H F 

Li metal mass balance oni Zee Re 

dr F 
LiCl mass balance antic = Rt- Rp, rz 1/kn , Ry = bro) — e amer, 

ot 1/kerkeo + t 
dT . 
Energy balance mcp T Aiset(Etn — E) — Aprh(T — Tamb) + Q 
Parameter Units Description 
(b) Symbols 

A m? Separator area (one side) 
AHT m? Heat-transfer area of cell 
D mie) Diffusion coefficient of LiAICl4 
c Equivalents m7? LiAlCl4 concentration 
Cp Tkg-!K7! Heat capacity of cell 
E V Cell potential 
En Vv Thermo-neutral potential for Li/SOCl 
F C eq ! Faraday’s constant 
h Wm? K! Faraday’s constant 
iset Am? Applied current per unit area of separator 
in x Am? Liquid-phase current density 
Jy Equivalents m7? SOCI; reduction rate per volume of cathode 
m Kg Mass of cell 
NL mol m~? Amount of lithium metal per area of separator 
NLicl mol m~? Amount of LiCl at negative per area of separator 
Q WwW Heat released by self-discharge reaction 
t s Time 
t None Anion transport number 
ty None Cation transport number 
T K Temperature 
Tamb K Ambient temperature 
Vo mi mol`! Partial molar volume of SOCI? 
Mio mi moll Partial molar volume of LiCl 
Vv mi mol7! Partial molar volume of LiAICly 
x m Distance from negative electrode surface 
a None Transfer coefficient 
E None Porosity 
di V Solid-phase potential 
di Vv Liquid-phase potential 
K Q27! cm7! Liquid-phase ionic conductivity 
v ms) Volume average velocity of electrolyte 
Ceff Q7 cm7! Effective electronic conductivity of cathode 
T None Tortuosity 
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ately after a discharge, the self-discharge was significant but 
decreased rapidly over time; on storage the self-discharge reac- 
tion decreased to a very low-rate, but could be reinvigorated by 
a pulse discharge. Nainville et al. [10] presented a percolation 
cluster model to describe passivation of a lithium anode that 
could predict the effect of additives on the thickness and poros- 
ity of the passive film. Here, a simpler model is presented for 
estimating the effect of the film on the discharge behavior with 
the goal of estimating lithium consumption over long periods of 
time. 


2. Model development 


Several authors have all presented very similar electrochem- 
ical models for the discharge of a thionyl chloride cell [2-4]. 
These models consist of equations describing mass balances on 
reactants, charge-transfer kinetics, porosity change in the cath- 
ode, electrolyte flow, temperature changes, and ohmic drops 
in the liquid and solid phases (see Table 1). These equations 
are explained in detail by Jain et al. [4] and their formula- 
tion was used in this work. This work extends previous models 
by accounting for direct reaction of the lithium electrode with 
thionyl chloride; that is, the self-discharge reaction. 

In the work mentioned above, Roth [9] found that the rate 
of heat generation from self-discharge, Q, had the functional 
form 


Q 


ao 


aay Á 


where t is time, ag and a; are constants, and x ranged from 1.17 

to 1.34. For consistency with Roth’s experimental results, the 

rate of formation of LiCl layer, R¢ (mol m~? s~!) was posited to 

decrease exponentially with the amount of LiCl, me (mol m~?) 
_ dnc 


Rr= = ke rue) 3 
f ep £0 (3) 


where be (mol m~? s~!)is arate constant and kr (m? mol!) isa 
constant. Integrating this equation and taking the time derivative 
gives 


= Lk 
(1/keikgo) + t 


Since the rate of heat generation due to self-discharge, Q (W) 
can be calculated as 


(4) 


A x A Hike 
Q=AxAH x R= L (5) 
(d/kfikfo) + t 


where A is active area (m7) of the lithium electrode. Eq. (5) is 
consistent with Roth’s Eq. (2) with x=1. 

Eq. (4) describes the formation of a LiCl film at the negative 
electrode. To account for Roth’s finding that pulse discharge 
accelerated the reaction between lithium and thionyl chloride, 
the following equation is posited to describe the destruction of 
the LiCl film. 


Ro = bat) — e Tune, (6) 


Table 2 
LR D-cell properties at build time 
Battery Property Value 
(a) Battery 
Voltage (V) 3.723 
Weight (g) 120.4 
Volume (cm?) 57.9 
Active area (cm?) 122.0 
Jellyroll property Value 
(b) Jellyroll 
Weight (g) 13.89 
Diameter (mm) 30.94 
Volume (cm?) 19.43 
Coating property Positive Negative 
(c) Active layers 
Thickness* (w/collector), (wm) 1803.4 1295.4 
Porosity (%) 83.5 0 
Total length (cm) 19.8 15.49 


è Overall thickness, since two-sided, thickness per active area is half. 


where Rp (mol m~? s~!) is the rate of film destruction, iy (A m~?) 
the current density, kpo (mol m2 s7!) a rate constant, and kp1 
(m? mol!) is a constant. This equation accounts for the destruc- 
tion of the LiCl layer by application of a current and the require- 
ment that LiCl must be present in order to be destroyed. 

Eqs. (4) and (6) constitute a kinetic model for LiCl film forma- 
tion and destruction at the lithium electrode. This kinetic model 
was incorporated into the electrochemical model described in 
detail by Jain et al. [4]. The equations constituting the electro- 
chemical model have been presented many times [2—5] and are 
summarized for reference here in Table 1. There is a slight differ- 
ence between the electrochemical model here and that used by 
Jain et al. [4] in that Jain used a complex empirical expression to 
represent the temperature dependence of the salt diffusion coef- 
ficient while in this work a simple Arrhenius equation was used. 
The partial differential equations defining the electrochemical 
model were solved using a pentadiagonal band solver. 

The code implementing the model was incorporated into 
Battery Design Studio® (see www.batdesign.com) in order to 
provide a user-friendly interface as well as provide capabilities 
for regression to experimental data and visualization of results. 


Table 3 

Discharge model parameters 

Parameter Values 
Positive thickness (single side) (um) 1220.59 
Do, salt diffusion coefficient (cm? s~!) 1.24e—02 
E,, diffusion coefficient (J mol~!) —3008.73 
ip, positive exchange current density (mA cm~?) 2.19e+06 
E,, positive activation energy (J mol~!) 5611.32 
io, negative exchange current density (mA cm~?) 2.65e+05 
E,, negative activation energy (Jmol~!) 4434.75 
Positive Bruggeman coefficient 2.26e—01 
Morphology factor 5.63e—02 
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Cell Voltage (V) 


Charge Removed (Ah) 


Fig. 1. Comparison of experimental [4] and simulated discharge curves at 
—18°C for loads 10, 50, and 250 Q. 


3. Results and discussion 


Table 2 lists the properties of the low-rate (LR) D-size cell. 
To verify the implementation of the model used here, the results 
presented by Jain et al. [4] were reproduced. Then the parameters 
determining the discharge behavior were refitted to the experi- 
mental data reported by Jain et al. [4] using the model described 
above (that is, simple Arrhenius expression for salt diffusion 
and self-discharge kinetics). In this work, a smaller value of the 
active electrode area was used (122 cm? versus 180 cm?) because 
a careful analysis of the cell construction revealed that only a 
fraction of the cathode area in the D-cell was actually facing 
lithium foil. Use of the smaller value necessitated refitting the 
parameters listed in Table 3. 

Table 3 summarizes the model parameters that best fit the 
experimental data. Note that the positive electrode thickness 
used to simulate discharge performance is significantly higher 
than the build value listed in Table 2. Although the electrochem- 
ical model does not account for the expansion of the cathode 
during discharge, the cathode is known to swell during dis- 
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Fig. 2. Comparison of experimental [4] and simulated discharge curves for a 
50 Q load at —55, —18, and 25°C. 
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Fig. 3. Comparison of experimental [4] and simulated discharge capacity as a 
function of the load over a temperature range of —55 to 49°C. 


Table 4 
Comparison of simulated mid-discharge voltages to experimental 
Discharge conditions Mid V Mid V Difference 
Resistance (Q) TCO) (simulation) (experimental) (%) 
50 25 3.57 3.54 0.99 
50 —18 3.29 3.30 0.39 
50 —55 3.01 2.97 1.10 
10 —18 3.10 3.20 3.12 
250 —18 3.47 3.43 1.26 
Table 5 
Comparison of simulated discharge capacities to experimental 
Discharge simulation Simulation Experimental Difference 
; (Ah) (Ah) (%) 
Resistance (Q) T (°C) 
50 25 14.98 15.6 3.92 
50 —18 10.33 9.97 3.69 
50 —55 3.38 3.52 3.84 
10 —18 4.27 4.44 3.63 
250 —18 14.99 15.4 2.64 
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Fig. 4. Comparison of simulated and experimental results [9] for results for 
self-discharge heat output for LR D-cell under open circuit conditions. 
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Fig. 5. Comparison of simulated and experimental results [9] for self-discharge heat output for LR D-cell after current pulses. 


charge [2,3,11]. Unfortunately, no data on the degree of swelling 
seems to be available to validate use of the fitted value. However, 
swelling of the cathode is necessary since the cathode pore vol- 
ume of the as-built cathode would limit cell capacity to 13.8 Ah 
(based on volume of solid LiCl formed) while experiments pro- 
duce >15 Ah discharge capacity. 

Figs. 1-3 compare simulation results to experimental data; 
Tables 4 and 5 give a quantitative comparison. The fit is compa- 
rable to that obtained by Jain et al. [4] and provides confidence 
that the model was correctly implemented. This model was then 
used to fit the calorimetric data of Roth [9]. 

For no-load, self-discharge conditions, parameters kpo and 
kp, are set to zero, and En and ke are adjusted; Fig. 4 compares 
experimental and simulation results. For pulse-load conditions 
that involve a series of single and continuous pulses of 1-s dura- 
tion, with duty cycles of 1%, 5%, and de (100%), kpo and kp, are 
adjusted; Fig. 5 compares experimental and simulation results. 
Parameter values giving the best fit are listed in Table 6. The 
simulation model gives a very reasonable fit to the behavior. 
This model was used to predict the effect of long-term (10 year) 
use of the cell. 

To assess the importance of lithium corrosion, an intermittent 
discharge was simulated where the rest period varied from 1 day 
to 30 days. A pulse current of 0.3 A of 5 s duration was applied 
and a rest current of 20 pA was used. Fig. 6 shows the con- 
sumption of lithium (in ampere hour) versus discharge capacity 
over a 10 year period. For the three cases considered, pulses 
every 1, 5, or 30 days, the total amount of lithium consumed by 


Table 6 

Kinetic parameters for Li/SOCly self-discharge (see Eqs. (3) and (5)) 
Parameter Value 

kro 5.21e—008 

key 2.8 

koo 2.95 

kot 3.6 


Ah Li Metal 


2.5 3 35 


0 0.5 l 1.5 2 
Discharge Capacity, Ah 


Fig. 6. Simulation 10 year intermittent discharge. Each line ends at 10 years. 


both self-discharge and discharge after 10 years is approximately 
8 Ah. So, there is no danger of depleting the lithium and thereby 
introducing a safety hazard. However, since the cathode-limited 
capacity is ~15 Ah, the lithium might be depleted if, at the end 
of 10 years, the cell were completely discharged. 

The simulation model presented here enables estimates of 
the discharge behavior of Li/SOCI, cells over a wide range of 
temperatures and discharge rates, as well as the self-discharge 
rate. This capability provides a rational basis for the design of 
battery-powered devices. However, the model is of limited use 
for design of batteries because it relies heavily on adjustable 
parameters (see Tables 3 and 6) to obtain good fidelity with 
experimental data. 
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